The Three- and Four-Nucleon Systems from Chiral Effective Field Theory 
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Recently developed chiral nucleon-nucleon (NN) forces at next-to-leading order (NLO) that describe 
NN phase shifts up to about 100 MeV fairly well have been applied to 3N and 4N systems. Faddeev- 
Yakubovsky equations have been solved rigorously. The chiral NLO forces depend on a momentum 
cut-off A lying between 540-600 MeV/c. The resulting 3N and 4N binding energies are in the same 
range as found using standard NN potentials. In additon, low-energy 3N scattering observables 
are very well reproduced like for standard NN forces. Surprisingly, the long standing ^4 H -puzzle is 
resolved at NLO. The cut-off dependence of the scattering observables is rather mild. 
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Effective field theory (EFT) has become a standard 
tool in modern physics and is applied to a large vari- 
ety of systems. It can also be used to construct nu- 
clear forces in a systematic and controlled manner. The 
spontaneously and explicitely broken chiral symmetry 
of QCD can be implemented in the EFT formulated in 
terms of the asymptotically observed Goldstone boson 
(pion) and matter (nucleon) fields. In the purely pi- 
onic and the pion- nucleon systems, there is an expan- 
sion parameter which is a typical external momentum 
(or the quark masses) divided by a hadronic mass scale 
of the order of the p meson or the nucleon mass. Any 
S-matrix element or transition current can be systemat- 
ically expanded in terms of this small parameter based 
on a systematic power counting. In systems with more 
than one nucleon, an additional non-perturbative resum- 
mation is mandatory to deal with the shallow nuclear 
bound states. This idea was put forward by Weinberg 
and taken up by van Kolck and collaborators in 
the construction of two- (NN) and three-nucleon (3N) 
forces. One basically constructs a potential based on 
the power counting and calculates bound and scatter- 
ing states by use of a properly regularized Lippmann- 
Schwinger or Schrodinger equation. One outstanding re- 
sult was that 3N forces (3NF) vanish to leading order |Q . 
Other groups also investigated low energy properties of 
NN systems along these lines A different counting 

scheme was proposed by Kaplan et al. Q (KSW) work- 
ing directly with the scattering amplitudes instead with 
an effective Hamiltonian like it is the case along the lines 
proposed by Weinberg. Another important feature which 
distinguishes the KSW approach from Weinberg's is the 
perturbative treatment of the one-pion exchange. Inde- 
pendent of these differences, such type of framework for 
the first time offers the possibility of calculating nuclear 
forces directly from fundamental principles and has a di- 
rect link to the chiral properties of QCD. Furthermore, 
this approach is firmly based on quantum field theory 
and avoids ill-defined concepts like meson-nucleon form 



factors. 

In 0] we have taken up Weinberg's idea and con- 
structed a NN and 3N potential based on the most gen- 
eral effective chiral pion-nucleon Lagrangian using the 
method of unitary transformations. In this method the 
field theoretical pion-nucleon Hamiltonian is decoupled 
such that an effective purely nucleonic Hamiltonian con- 
sistent with a power counting scheme arises. Previous re- 
sults were obtained in time-ordered perturbation theory, 
which lead to energy-dependent nuclear forces. In our 
formalism, we arrive at hermitian energy-independent 
nuclear forces which we consider to be a major advantage 
with respect to applications to 3N and 4N systems, the 
issue of this letter. In |l|,||j7[|^] NN forces have been de- 
veloped at leading, next-to-leading and next-to-next-to- 
leading orders, LO, NLO and NNLO, respectively. At LO 
the potential is represented by the ordinary one-pion ex- 
change (with a point-like coupling) as well as two contact 
interactions without derivatives. At NLO one includes 
the leading chiral two-pion exchange as well as all possi- 
ble contact interactions with two derivatives, whereas at 
NNLO we have additional two-pion exchange with low- 
energy constants (LECs) determined from pion-nucleon 
scattering The forces are properly renormalized and 
contain nine parameters related to those four-fermion 
contact terms. The one- and two-pion exchange pieces 
are parameter free. The nine LECs have been uniquely 
fixed to low energy NN phase shifts in the s- and p- waves. 
The parameter free predictions for higher energies and 
partial waves and also deuteron properties are in gen- 
eral rather good. It was also observed that the NNLO 
predictions are better than the ones based on the NLO 
potential, as expected in a systematic EFT. 

The natural question arises now, whether the NN 
forces based on chiral perturbation theory (%PT) will be 
also successful in describing 3N and 4N low-energy ob- 
servables. To that aim we solve the Faddeev-Yakubovsky 
equations rigorously for 3N and 4N systems pfj|-Jl2[ and 
determine binding energies and various scattering observ- 
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TABLE I. Neutron-proton phase shifts in our approach 
(upper row) compared to the Nijmegen PSA (middle row) 
and the CD-Bonn potential (lower row). 
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ables. To the best of our knowledge this is the first time 
that xPT has been practically applied to nuclear systems 
beyond A — 2 within the Hamiltonian approach. 

In this first application we restrict ourselves to the 
NLO NN potential. In a forthcoming article we shall 
go on to NNLO and include also 3NFs, which occur at 
that order the first time. The NLO results presented here 
are therefore parameter free and can serve as a good test- 
ing ground for the usefulness of the approach. Of course, 
some aspects of the 3N system have already been stud- 
ied in nuclear EFT (l^Jl^], but not as direct extensions 
of the NN system as done here. 

In order to use the chiral NN forces in the NN 
Lippmann-Schwinger equation one has to introduce a 
momentum regulator A. We remark that this regulariza- 
tion on the level of the scattering equation is completely 
different from standard methods which are applied to in- 
dividual diagrams. Here we use the smooth regulator, its 
precise form is given in [^J. In order to investigate the 
cut-off dependence of 3N and 4N observables we have 
generated several NN potentials corresponding to differ- 
ent cut-offs between A =540 and 600 MeV/c. They were 
all fitted to the l S , 3 S 1 - 3 D 1 , 1 P 1 and 3 P a,2 NN phase 
shifts up to £i a b = 100 MeV (for the potential parame- 
ters contact E.E.). In || we had already demonstrated 
that going to higher order in the EFT reduces the cut- 



TABLE II. Theoretical 3 H and 4 He binding energies for 
different cut-offs A compared to CD-Bonn predictions and to 
the experimental 3 H binding energy and the Coulomb cor- 
rected 4 He binding energy. The kinetic energies and S' , P 
and D state probabilities for 4 He are also shown. 
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off dependence and allows to choose larger values for the 
cut-offs, as expected from general arguments [ p~5| . The 
resulting phase shifts for NLO are compared to the Ni- 
jmegen phase shift analysis p6| and the ones obtained 
from the CD-Bonn potential fll7| ] in Table |. The agree- 
ment is fairly well and we know from Q that one has to 
go to NNLO to improve systematically on that. Also, we 
restrict ourselves to the isospin symmetric case. Thus we 
do not take into account various charge independence and 
charge symmetry breaking effects like e.g. the pion mass 
difference, which are known to be significant at very low 
energies. Such effects can also be included systematically 
in our EFT @. 

Let us regard 3N and 4N binding energies now. In ad- 
dition to the relative motion in the NN subsystem there 
occur relative motions of the third and fourth nucleon 
now. In the spirit of a low momentum theory one has 
to expect that those additional relative momenta should 
be also small. This turns out to be true as the numerics 
tells us. For standard NN potentials, which are defined 
for all momenta up to infinity, one introduces a momen- 
tum cut-off for the numerical integrations. Since the 
pertinent integrals are convergent, the results are cut-off 
independent. Using the chiral NN potential we could re- 
duce the cut-offs in those relative momenta from their 
typical values 8 fm _1 arising with standard potentials 
down to 4 fm -1 without changing the result. In the 
Faddeev-Yakubovsky equations there is no mechanism 
which could introduce high momenta if the NN subsys- 
tem momenta are small to start with. We find for the 
fully converged solutions the 3N and 4N binding energies 
as given m Table g. The ranges are compatible with 
what is found using realistic potentials ju^. Note that 
the NN forces were included up to the total NN angular 
momentum of j max =6. 

Also 3N scattering can be solved rigorously nowadays 
|]l9| and we show in Figs. 0, || a small selection out of 
the great wealth of observables in comparison to data 
and the theoretical predictions of CD-Bonn. Three ener- 
gies, one below the nd break-up threshold and two above 
are chosen. Note further that we do not indicate error 
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FIG. 1. (a): The nd- total cross section (in barn) for the 
chiral forces (A = 540 MeV/c, dotted curve; A = 600 MeV/c, 
dashed curve), and CD-Bonn (thick solid curve), (b), (c) 
and (d): Chiral NN force and CD-Bonn predictions for nd 
break-up cross section da = dn f d ^ dS [^jtv] at £i ab =13 
MeV along the kinematical locus S. The various break-up con- 
figurations are chosen according to Figs. 42, 39 and 35 in [ 19 1, 
respectively, pd data are (v) pjl; nd data are (a) 0], (•) [ 22 1 , 
(♦) & 



bars for the data since in most cases they will be not 
distinguishable on this scale. In all cases we show the 
predictions of the chiral NN potentials for the cut-offs 
A = 540 and A = 600 MeV/c and compare it to the pre- 
diction of one of the most modern, so-called realistic NN 
potentials, the CD-Bonn p7| . As the simplest observ- 
able we show first the nd total cross section in Fig. |l](a). 
The three theoretical curves overlap completely at very 
low energies and then the chiral predictions start to de- 
viate somewhat from the data as expected for our EFT 
at NLO. 

In case of the 3N break-up we selected in Fig. |l](b,c,d) 
a few often investigated configurations, the space-star, 
a final state interaction peak configuration and a quasi- 
free scattering (QFS) configuration, respectively. We find 
very good agreement of the chiral NN force predictions 
with the one from CD-Bonn. In case of QFS (d) and 
space-star (b) some of the discrepancies are expected to 
be caused by Coulomb force effects not included in the 
theory. The upper group of data in (b) are nd data and 
the disagreement with the theory presents a well known 
puzzle at the moment |p4|j . 

For elastic nd scattering we display in Fig. || the an- 
gular distributions, the nucleon analyzing power A y and 
the tensor analyzing powers Ti^. Except for A y there 
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FIG. 2. nd elastic scattering observables at E n — 3MeV 
(left column) and E n = lOMeV (right column), pd data are 
Mi @, (v) @, (a) @, (o) @. nd data are (a) [f§, (▼) 
[p0|. For remaining notations see Fig. hi. 



are no nd data for those energies. The discrepancies be- 
tween data and theory for T-i^ and for the differential 
cross section can be traced back to pp Coulomb force ef- 
fects ||l| . Thus except for A y the agreement of CD-Bonn 
(thick solid curve) with the data is rather good, which is 
a well known fact and is just given for the sake of ori- 
entation. The dotted and dashed curves are the chiral 
NN force predictions for A = 540 and A = 600 MeV/c, 
respectively. The da and Tik agree rather well with the 
CD-Bonn result and thus with the data. We consider this 
to be an important result, demonstrating that the chiral 
NN forces are very well suited to also describe low-energy 
3N scattering observables rather quantitatively. On top 
of that, surprisingly for us, the chiral force predictions 
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can be used directly instead of wavefunctions generated 
by standard NN forces. This appears advisable since in 
this study for NLO we noticed that kinetic energies in 3N 
and 4N bound states are smaller and two-nucleon correla- 
tion functions are smoother than for standard NN forces, 
indicating a change of the wave functions. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft (A.N.) and by the Polish Committee for Sci- 
" 60 120 180 " 60 120 180 entific Research under Grant No. 2P03B02818 (H.W.). 

6 [cleg] e [cieg] The numerical calculations have been performed on the 

FIG. 3. np analyzing powers A y at 3 (left panel) and 12 CRAY T90 and T3E of the John von Neumann Institute 
MeV (right panel) for the chiral forces (A = 540 MeV/c, for Computing in Julich, Germany, 
dashed curve; A = 600 MeV/c, dot-dashed curve), PSA (solid 

curve) and CD-Bonn (dotted curve). 




are now significantly higher in the maxima of A y than 
for CD-Bonn and break the long standing situation that 
all standard realistic NN forces up to now underpredict 
the maxima by about 30 %. This is called the A y -puzzlc 
p9[ . We are now in fact rather close to the experimental 
nd values. Since we restrict ourselves to NLO we can not 
expect a final answer from the point of view of chiral dy- 
namics but this result for A y is very interesting. In that 
context it is important to note that on a 2N level the chi- 
ral potential predictions agree well with the predictions 
based in the Nijmegen PSA as is shown in Fig. ||[ There 
we also include the CD-Bonn predictions. This agree- 
ment is especially important for the np analyzing power 
A y , since it is rather sensitive to the 3 Pj NN phase shifts, 
which influence also strongly the 3N A y . 

These very first results using chiral NN forces in 3N 
and 4N systems are very promising. Since we restricted 
ourselves to NLO we could not expect a very good de- 
scription of the data, since on-shcll properties are not 
perfectly well fitted. But the results show that these ef- 
fective chiral forces are very well suited to describe also 
low energy properties of nuclear systems beyond A = 2. 
They agree rather well with standard nuclear force pre- 
dictions as exemplified with CD-Bonn and most impor- 
tantly they break the stagnation in the A y puzzle. Our 
result provides a counter example to the suggestion |3^] 
that NN forces alone can not describe 2N and 3N A y 's at 
the same time and 3NF's should cure the 3N A y puzzle. 
Examples for such trials can be found in pl|p2[ . 

It will be very interesting to perform the next step and 
use the NNLO NN forces, which is a systematic improve- 
ment. On this level also 3N forces have to be incorpo- 
rated the first time, which in xPT are defined consis- 
tently with the NN force. It should be mentioned further 
that due to the underlying Lagrangian the coupling of ex- 
ternal probes (photons for instance) is well defined and 
exchange currents consistent to the nuclear forces can be 
generated. Also the hybrid approaches |m|,[35) can be 
avoided and wavefunctions based on the chiral dynamics 
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